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Abstract 

Infinitely many explicit solutions of certain second order differential equations with 
^ ■ an apparent singularity of characteristic exponent —2 are constructed by adjusting the 

parameter of the multi-indexed Laguerre polynomials. 
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In recent years, global solutions of Fuchsian differential equations with more than three 
regular singularities and those of the differential equations of confluent types, attracted 
attention [T] . The multi-indexed Laguerre and Jacobi polynomials [21 [3] and their simplest 
(one-indexed) cases, i.e. the exceptional orthogonal polynomials, provide explicit examples 
of infinitely many global solutions of such equations. They are obtained by Darboux-Crum 
transformations [I] [6] applied to the exactly solvable Schrodinger equations [8] of the 
radial oscillator potential and the Poschl- Teller potential. In these global solutions, all the 
extra singularities turn out to be apparent and the characteristic exponents are all — 1. 
This is a generic property of the Darboux-Crum transformations. In the present paper, we 
construct, by adjusting the parameter of the multi-indexed Laguerre polynomials, several 
families of infinitely many global solutions of Schrodinger equations with two or more finite 
regular singularities and one irregular singularity at infinity. At one of the regular singular 
points, these solutions have an apparent singularity of characteristic exponent —2. We now 



show one of the examples which is labeled as (E) in section three. The Hamiltonian is given 
by 

u _ * , x 2 , 195 16 2 ^log(x 2 + 6) 3 (x 2 + 14) ^ 
dx 2 Ax 2 dx 2 

We change the variable by 77 = x 2 . Then the equation 1-Lip{x) = Sip(x) (S G C) in the 
variable 77 has apparent singularities rj = —6 and 77 = — 14, whose characteristic exponent 
is —2 and —1 respectively. All the apparent singularities are located outside of the half 
line (0,oo), thus the eigenf unctions of the Schrodinger equations are constructed by using 
Darboux-Crum transformation, and they are global solutions of the differential equations. 
The essential part of the numerators of these global solutions form orthogonal polynomials 
over a half line (0, 00), with the square of the denominator forming the part of the weight 
functions e~ x x 7 /{(x + 6) 4 (x + 14) 2 }. The polynomials in the numerators of the global 
solutions are obtained within the framework of multi-indexed orthogonal polynomials, and 
their degrees have "gaps". In this example, the polynomials of degrees 0,1, and 2 are missing. 
We also report examples (A), (B), (C), (D), (F), (G) in the two index case and an example 
in the three index case. Each example has its peculiar properties (see section three). The 
orthogonal polynomials related to (A) and (B) are also reported in [9]. 

This paper is organised as follows. In section two, we review the Darboux transformations 
with the induced singularities. In section three we show that the special cases of multi- 
indexed Laguerre polynomials provide several families of infinitely many global solutions of 
confluent differential equations with two or more finite regular singularities and one irregular 
singularity at infinity. They have an apparent singularity of characteristic exponent —2. The 
final section is for a summary and comments. 

2 Darboux transformations and induced singularities 

Our starting point is a generic Schrodinger equation with a Hamiltonian % 

n = ~ + u(x), (2.1) 

in which U(x) is a meromorphic potential. Let ip(x) and <p(x) be two distinct solutions, not 
necessarily square-integrable eigenfunctions, of the Schrodinger equation: 

U^{x)=£^{x), H<p(x) = £tp(x), £,£eC. (2.2) 
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Then it is elementary to show that a new function 

^ f W ^^]( x ) _ v( x W( x ) ~ <P'( x )i>( x ) 
<p(x) <p(x) 

is a solution of a new Schrodinger equation with a deformed Hamiltonian 71^ 

= + ^g^j.jf^W, (2 .4) 

with the same energy £: 

H^ij {1) (x) =S^ 1 \x). (2.5) 

The zeros of the "seed" solution ip(x) produce singularities of the equation, and there are no 
multiplicities of the zeros of <p(x) except for the singularities of U(x), because the function 
(p(x) satisfies the Schrodinger equation and the exponents at a holomorphic point is or 1. 

To obtain the deformed Hamiltonian which admits an apparent singularity of exponent 
—2 or higher, we use multiple Darboux-Crum transformations [I] [6]. Let {(fj(x), £j}, j = 
1, . . . , M be distinct seed solutions of the original Schrodinger equation (12.21) and also ip(x) 
be a solution: 

Hil>(x) = £i/j(x), Wpj{x) = £j<fj(x), £, £j e C, j = 1, . . . , M. (2.6) 
By repeating the above Darboux transformation M-times, we obtain a new function 

V%i,...,y3 M ](x) 

which satisfies an M-th deformed Schrodinger equation with the same energy [5],[6]: 

U (M W M \x) = £^ M \x), (2.8) 

H M = * +U M {X)} ^( g )g^)-2 jil0gWly ; i ;- ytfl(z) , (2.9) 
ax 2 ax 2 

M 

(ij( M \^) = l[(£-£ J ).^^). (2.10) 
Here W[/i, . . . , /jv](x) is a Wronskian 

/ l<j,k<N 

and the inner product among real functions is defined as usual: 

if,g) = / f{x)g{x)dx. 



Except for the singularities of the potential, a zero of the Wronskian W[v?i, . . . , <Pm](x) 

W[<p u <p M ](x) = c (x - x ) m + 0(x - x ) m+1 , m G Z >0 , (2.11) 
corresponds to an apparent singular point of the deformed Hamiltonian W M ^: 

x ps Xq, U^ m Mx) = U(xo) + 7 + regular terms. (2.12) 

(X-Xo) 2 

The characteristic exponents of the solutions of the new Schrodinger equation (12. 5 j) at the 
apparent singularity x are determined by 

p{p-l)-2m = 0. (2.13) 
For certain values of the integer m, the monodromy of the added singularity becomes trivial: 

m=l, p = -l, 2, (2.14) 
m = 3, p=-2, 3, (2.15) 
77i = 6, p = -3, 4, etc. (2.16) 

When all the added singularities are of trivial monodromy, the possibility of global solutions 
arises. The infinite families of the exceptional (or one-indexed) Laguerre and Jacobi polyno- 
mials [2] are the examples of the global solutions of second order Fuchsian (and its confluent 
limit) equations with many extra singularities with all p = — 1. This is the generic case, since 
any solution (p(x) of the Schrodinger equation (12. 2p cannot have a multiple zero (m > 2). 

3 Multi-indexed Laguerre polynomials 

In this section, we present several families of infinitely many global solutions with a high 
degree (p = —2) of apparent singularity. The strategy is simple. We start from certain 
"non-generic" multi-indexed Laguerre polynomials and choose the parameter g such that a 
triple zero is achieved for the seed solution (the Wronskian). Then we construct the global 
solutions explicitly. When the zeros are outside the half line (0, oo), the global solutions are 
square integrable and they are eigenfunctions of the deformed Hamiltonian. 

3.1 "Non-generic" multi-indexed Laguerre polynomials 

For the multi-indexed Laguerre polynomials [2], we adopt the Schrodinger equation (12.21) 
with the radial oscillator potential [8] 

U(x; g) = x 2 + 9 ^ 9 ~ ^ - (1 + 2g), g > 0, < x < oo. (3.1) 
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It has a regular singularity at x = with the characteristic exponents p = g, 1 — g and an 
irregular singularity at x = oo. The eigenpolynomial solutions are 

"H0 n (x; = £ n 0„(x; g), £ n = An, (3.2) 

^.(il^e^Vli" 5 '^)), V (x)=x 2 , n = 0,l,..., (3.3) 

in which Lffl(rj) is the degree n Laguerre polynomial in 17. They form a complete set of 
orthogonal polynomials over (0, 00): 

d f f°° 1 /"°° (9 ~) (9 ~) 

(0n,0m) = / 4> n (x; g)(f) m (x; g)dx = - / e' r, r] 9 ^Ln 2 (77) Lm 2 (r/)rf7/ (3.4) 
= h n (g)5 nm , h n (g) = —T(n + g + \). (3.5) 

The generic multi-indexed Laguerre polynomials are constructed by the multiple Darboux- 
Crum transformations (12.71) in terms of two types of polynomial (v e Z>o) seed solutions, 
called the virtual state solutions [2], 

^(^^^e^VL^^^x)), Sl(g)^-A(g + y+l), (ft, ft = 00, (3.6) 
^'(^l^e-^x^Wx)), ^^^-v-i), (ft, ft) = 00, (3.7) 

with certain constraints on the range of the parameter g. In order to generate a wider class 
of global solutions with a high degree of apparent singularities, we will adopt non-virtual 
state polynomial seed solutions, too: 

^{x;g)*e&x l -'I$~'\-v(x)), ^ n (s) d = f -4(v+l), (ft 1 , ft 1 ) = 00. (3.8) 

Here we summarise the definition of "non-generic" multi-indexed Laguerre polynomials. 
Definition: The "non-generic" multi-indexed Laguerre polynomials are the main part of 
the eigenfunctions ^ M \x) (12. 7p generated by the multiple Darboux transformations with 
the eigenfunctions {</> n (x; g)} and the three types of seed solutions (I3.6p . ( 13 .7p and (13. 8p . the 
virtual and non-virtual state solutions, and that the parameter g is restricted only by the 
positivity g > 0. 

The regularity of the solutions and/or the positivity of the resultant weight functions 
etc. must be verified in each specific case for the chosen parameter values. This is in good 
contrast with the multi-indexed Laguerre and Jacobi polynomials defined in [2] , in which the 
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parameters are restricted (eqs.(23) and (24) of [2]) so that the seed solutions are of definite 
sign and their inverse are not square integrable 



(l/#, l/#),= (l/0j,l/^)=oo. 

These two additional properties guarantee the regularity of the solutions and the positivity of 
the weight function. In other words, these restrictions are lifted for the possible construction 
of the non-generic multi-indexed Laguerre polynomials. In the singular solution cases, the 
square integrability is irrelevant and we also consider the g < example (F). 

For the M-index case, the seed solutions is now W[</? Vl , . . . ,(p VM ](x). We look for the 
possibility that this function has a triple zero or higher zeros at some points. 



3.2 Two- index case 



Here we consider the simplest multi-index case, the two-index case: 



W[y?i, ip 2 ,<Pn](x) 



(3.9) 



We focus on the seed solutions of the form 



W[4>^(x;g),4>l 2 (x;g)](x), W$™(x;g),$ 2 (x;g)](x), W[^ vi (x;g),^ 2 (x;g)}(x), 

with lower degrees v x and v 2 and evaluate the discriminants of the polynomial part of the 
seed solutions. At the roots of the discriminant, which is a polynomial in g, the above seed 
solution has a multiple zero. The first two are the non-generic type and the last one is the 
generic type. As we will show in some detail, they have different features. 
So far we have encountered the following seven cases: 

(A): W[ff(s;f), !)](*), 



(C): W[^(x;}),^?(x;5)](x), 
(E): W[^(x;f),^(x;f)](x), 



(B): W[$\x-,\),fc(x-,\)](x), 
(£>): W[ff(s;S),®s; })](*), 
(F): W[4>\(x;-f),^(x;-f)](x), 



There could be some others which we have missed. 



6 



3.2.1 Seed solution case (A) 

In this case 

W[4F(x;g),4>l(x;g)](x) 

= -^(2g + 1) e x2 (-9 + 18g + Ag 2 - 8g 3 + (18 - 8g 2 )x 2 + (12 + 8g)x 4 + 8x 6 ) . 
The discriminant of the polynomial in the variable r](= x 2 ) is 

2048(2^ -3)(22 + 3) 2 (%-3) 2 , 

and g — | gives a cubic zero 

W[ff(x;f),S(x;|)](x) = ^e* 2 (3 + 4x 2 ) 3 . (3.10) 

The potential is 

„ . , 2 3 48 288 13 

U A (x)=x 2 - + 5 . 3.11 

V ; 16x 2 3 + 4x 2 (3 + 4x 2 ) 2 2 V ; 

3.2.2 Seed solution case (B) 

In this case 

W{$\x;g),fc(x;g)}(x) 

= -^(2g - 3) e x2 (5-2g- 20g 2 + 8g 3 + (10 + 16g - 8g 2 )x 2 + (20 - 8g)x 4 + 8x 6 ) . 
The discriminant of the polynomial in rj is 

-2048(2( ? -5) 2 (2( 7 + l)(4( ? -l) 2 , 

and g — \ gives a cubic zero 

W[ff(z;±),#(s;±)](a;) = (3 + 4a; 2 ) 3 . (3.12) 

The potential is the same as Ua(x) up to an additive constant: 

rr / X rr / N 2 3 48 288 11 . 

!/,(,) = + W - _ + _ - ^—^ - — . (3.13) 
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3.2.3 Seed solution case (C) 

In this case 
W[ffi(x;g),$(x;g)](x) 

= \x 3 ~ 2a (-135 + 17Ag - 68g 2 + 8g 3 + (-54 + 4% - 8g 2 )x 2 + (36 - 8g)x 4 + 8x 6 ) . (3.14) 
8 

The discriminant of the polynomial in 77 is 

-2048(2(7 - 9)2(2(7 -3)(4^-9) 2 , 

and g — | gives a cubic zero 

2\3 



W[^(x;S),0f(x;S)](x) = v ,L,/ 2 • ( 3 - 15 ) 



(3 + Ax 
6Ax 3 / 2 

The potential is exactly the same as £/#: 

rT , ^ tt , ^ 2 3 48 288 11 

Uc(x) = Ub\x) = x H 77 . 

CK J K ' 16x 2 3 + Ax 2 (3 + 4x 2 ) 2 2 

Since the eigenvalue £ n = An is independent of g, the corresponding solutions </>n (x) ( 13. 9ft 
for (B) and (C) satisfy the same equations. Thus they are identical. We will not consider 
the case (C) in the next subsection. 

3.2.4 Seed solution case (D) 

In this case 

W[ffi(x;g),$(x;g)](x) 

= - x 3 - 2g (-135 + 17Ag - 68g 2 + 8g 3 + (54 - A8g + 8g 2 )x 2 + (36 - 8g)x A - 8x 6 ) . 
8 

The polynomial part is obtained from that of (I3.14p by the change x 2 — > —x 2 . Thus the 
discriminant is the same and g = | gives a cubic zero 



[— % 4- At 2 ) 3 

W[^\x; J),^(x; l)](x) = 6 ^JJ . (3.16) 



The potential has a singularity in (0, 00): 

rr / X 2 3 48 288 11 , „. 

= X ~T^ + 33Ttf + (-3 + ~ T (3 - 17) 

The corresponding solutions (13. 9p are not square integrable and they do not form orthogonal 
polynomials. 

The cases (E)-(G) belong to the generic type. 
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3.2.5 Seed solution case (E) 

In this case 

W[4>l(x;g),^(x;g)](x) 



The discriminant of the polynomial in r] is 

-16(20 - 15) 2 (2£ - 3) (2<? + 1)(2<? + 3) 2 
and g = y & Yes a cubic zero 

Wfe f ), f )](*) = (6 + * 2 ) 3 (14 + x 2 ). 

The potential is 

2 195 144 12 112 4 

3.2.6 Seed solution case (F) 

In this case 

W[4>\(x;g),$(x;g)](x) 

= — -g- — + 4g 3 - g 4 - (— + 7g- Ug 2 + 4g 3 ) x 2 



16 2 V 2 

- ( - - 160 + 6g 2 j + (6 - Ag)x 6 - x 8 . 



2 

The discriminant of the polynomial in rj is 

-16(20 - 5) 2 (20 - 3) (20 + 1)(20 + 13) 2 

and at a negative value of the parameter g = — y it gives a cubic zero 

W[#(s; -?), -f )](*) = -(* 2 - 14)(:r 2 - 6) 3 . 

The potential is singular in (0, oo): 

TT . . 2 195 144 12 112 4 



3.2.7 Seed solution case (G) 

In this case 



W[%(x;g),ffi(x;g)](x) 

= 1 {(-3 + 2g)(-l + 2g)(l + 2«?)(3 + 2g){5 + 2g) 

+ 10(-3 + 2g){\ + 2#)(3 + 2g)(h + 2g)x 2 + 80(-l + 2g)(3 + 2g)(5 + 2g)x i 
+ 160^(5 + 2g)x 6 + 80(3 + 2g)x 8 } . 

The discriminant in r\ is proportional to 

(2g - 3)(2g + l)(2g + 3) 2 (2g + 5f(10g - 39) 2 

and g = || gives a cubic zero: 

Wfef ),0r>;| )](x) = ^(12 + 5x 2 ) 3 (2244 + 495x 2 + 25x 4 ). (3.22) 

The potential is 

1131 44 1440 60 



U G (x) =x 2 + 



100x 2 5 (12 + 5a; 2 ) 2 12 + 5a; 2 

165000a; 2 20 ( lO.r 2 - !)!) ) 

+ (2244 + 495a; 2 + 25a; 4 ) 2 + 2244 + 495a: 2 + 25a; 4 ' 



The basic structure of the (G) case is essentially the same as that of (E). So we will not 
present the explicit formulas of the (G) case, since they tend to be rather lengthy. 

3.2.8 The underlying structure of these seed solutions 

In the non-generic cases (A) and (B), the potential functions are essentially the same, but 
the main eigenf unctions carry different coupling constants, {(f) n (x; |)} and {(p n (x; |)}. Thus 
the corresponding multi-indexed polynomials are entirely different. The equality (up to a 
sign) of the seed solutions of (A) and (B) is easily understood, when one notes the trivial 
identity 

4>l(x;g)=4>™(x;l-g). 

In fact, we have 

\) = 0i n (x; |) = ^ 2 a;i(3 + 4x 2 ), (3.24) 
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(1/fe J), Vfe J)) = (l/fffo |), l/0i n (x; f )) < oo, (3.25) 
fe f) = i) = ^e^/ 2 a;I(3 + 4a; 2 )(15 + 4* 2 ), (3.26) 

|), l/£(x; f )) = (l/fffo J), l/fffa J)) = oo. 

A cubic zero of the seed solution in (A) and (B) is realised by a nice magic of the Wronskian. 
The fact that l/(f>\(x] |) and l/^ l (x; |) are square integrable f!3.25j) means the departure 
from the mechanism of multi- index polynomials developed in [2] . 
For the cases (C) and (D), the functions are: 

*?(*;!) = -^ 2 / 2 x-i(3 + 4x 2 ), = ^e-^ 2 x-i(-3 + 4x 2 )(l + 4x 2 ), (3.27) 

fffol) = ^' 2 / 2 x-I(-3 + 4x 2 ), ff(x; J) = ^ 2 / 2 x-I(3 + 4x 2 )(-l + 4x 2 ), (3.28) 

|), |)) = oo, (l/fffo f ), l/0?(a;; |)) = oo. 

Again the cubic zeros of these cases are realised by a nice magic of the Wronskian. As 
expected, <$(x; |), 0i n (x; |) and 2 n ( x > |) change sign in (0, oo), which is a deviation from 
the mechanism of multi- index polynomials developed in [2]. 

Now it is clear that the solutions (fr n (x) (13. 9 p have a p = — 2 apparent singularity, that 

is, 

« (3+ 1 4x2)2 , (A), (5), (C7), 0i 2) (x) cc F 3^ I , (U). (3.29) 

Among the six terms of the numerator Wronskian W[(pi,ip2,<j ) n\{x) of (13. 9p . the four terms 
containing <j>' n {x) or (p'^{x) contain at least one term (p±(x) or (p2(x) without a derivative. 
Thus they contain at least one factor of 3 + 4x 2 in (A), (B), (C) and —3 + 4x 2 in (D). The 
remaining two terms which contain <p n (x) is proportional to 

(p[(x)(p%(x) - <Pi(x)(p' 2 (x). 

Since (pj(x) is a solution of the original Schrodinger equation ( 12. 2D . they have 

<p';(x)=(u(x)-£j)<pj(x), j = 1,2. 

Thus the above terms also contain one factor of <p\{x) or (p2(x), that is, 3 + 4x 2 in (A), (B), 
(C) and —3 + 4a; 2 in (D). Thus one factor of the cubic zero of the seed function W[<pi, y^K^) 
is cancelled and we arrive at (13.291) . 
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The relationship between the generic cases (E) and (F) has similarities with that between 
the non-generic cases (C) and (D). The functions in (E) are 

0i( x; f) = I e x2/2 xT(6 + x 2 )(12 + x 2 ), ${x; f ) = -e^^f (6 + x 2 ). (3.30) 

The functions in (F) are 

-f ) = e a;2/2 a;^T(-6 + x 2 ), -f ) = ^e~ x2/2 xT(-6 + x 2 )(-12 + x 2 ). (3.31) 

Thus it is also clear that the solutions (fy n (x) ( \3.9\j have a p = — 2 apparent singularity, that 
is, 

tf)(s)oc (6T^' K (3 ' 32) 

3.3 Global Solutions with a p = —2 apparent singularity 

Here we present the infinitely many global solutions with a p = —2 apparent singularity 
in explicit forms. They are simply obtained by evaluating (jr n (x) ( 13. 9p in each case ( A)- 
(F), to be denoted by 4> A (x), <f>^(x),. . . , 0^(x), respectively. Most formulas look simpler as 
functions of rj = 77 (x) = x 2 , since the eigenpolynomials of the radial oscillator potential are 
the Laguerre polynomials in rj, ( 13. 3ft . 

3.3.1 Global solutions (A) 

The global solutions are 

-H A ^(x)=S n ^(x), n A = -^ + U A (x), S n = 4n, (3.33) 

ax z 

-77/2 3/8 

tf( g )= + » = °> ( 3 - 34 ) 
in which C A {rj) (n > 0) is a degree n + 3 polynomial in 77 defined by 

££(77) d ^ f (-117 + 15677 + 20877 2 + 6V) 4 1/4) fa) 

- 4n(3 + 4r7) 2 L( t 1/4) (77) - 477(3 + 4^)(15 + 477) (77). (3.35) 

In evaluating the Wronskian W[0* n (x; f ) , ^(^j |))^n(^j f)]^); the second derivative of the 
Laguerre polynomial Ln (tj) appears. It is replaced by Ln^\rj) and d n Ln (tj) by using 
the equation for the Laguerre polynomial: 

xd 2 x L^ (x) + (a + 1 - x)8,4 Q) (s) + «4 q) (s) = 0. (3.36) 
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There is another member of this family corresponding to the eigenvalue S = —8, be- 
cause of the fact that the Darboux transformation in terms of 0i H (x; |) is not one to one, 
(l/^tx; |)> Vffta !)) < fl3T25|) . see 0. Let us denote it tentatively by C A 2 (r]): 

£^^15 + 477, (3.37) 

whose normalisation is irrelevant. Note that n = — 1 is missing. Some lower members are 

765 

££(77) = -117 + 15677 + 208r? 2 + 64r? 3 , Cf(r]) = — — + 408r? + 408r/ 2 - 64r? 4 , 

Ci(n) = + ^77 + 5Q7 V 2 - 252?7 3 - 168r/ 4 + 32r/ 5 , 

32 8 

„ A . , 43875 20475 4875 2 ^ rn 3 orn 4 400 5 32 6 

£f(v) = -^28" + "le - ^ + "T" 77 ~ v ~ v + T v ~ Y v ' ■ " ■ " 

They are also discussed in §6.2.5 of [5]. They satisfy the orthogonality relation 

1 

, 3+ ' Ci(r,)C&(v)dr, = n, m = -2, 0, 1, . . . (3.38) 

g 4(n + 2)(4» + 18)r(n + i) ^ _ 2 
n! 

which is obtained by rewriting (I2.10p with £{ n (§) = —8 and ^(f) = —13. 
By rewriting the Schrodinger equation for 4>^(x) (13.331) in terms of 

3 

<jfi{x) = e-" /2 r/8 7/(77), 

we obtain a second order differential equation with regular singularity at 77 = and 77 = — | 
and an irregular singularity at x = 00: 

W" + (~ ~ V)V' + 45 "3 ^^y^ y + ny = 0, n = -2,0,l,...,. (3.40) 
For each n we have a global solution with a p = —2 apparent singularity: 



(3 + 4T7) 2 



Vn(v) = -o , A % > n = -2,0,l,... > . 



The differential equations for the polynomials {£„(??)} read 



477(3 + 477)^(77)" - (477 - 1)(15 + 477)^(77)' 

+ 4(5 + 12r/ + 77(3 + 4t7))£^(?7) = 0, 77 = -2, 0, 1, . . . , . (3.41) 
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The groundstate wavefunction 0^ 2 ( x ) Q3.34j) with C A 2 (13.371) and the energy £ = — 8 
is related to the potential U A (x) (13. lip through the prepotential w A (x) (Riccati form of 
Schrodinger equation): 

W») def e-" 2 / 2 xf (15 + 4x 2 ) f dw A (x) \ 2 d 2 w A (x) 

6 = — (sT^y — * + ^x^~ 8 = Ua{x) - 

In other words, the Hamiltonian H A (13.331) is written in a factorised form: 

H A = A«A A -%, A A^d_dwA^^ AA] = _d__d_w^_ ^ 

doc doc d/Oc d/Oc 

It is easy to verify that the standard Darboux-Crum transformation [5] of deleting the 
groundstate applied to (I3.42[) will lead to the system of one-indexed (or the exceptional) 
polynomial system generated in terms of the virtual state solution (p\(x; |): 

H M def jAjA] g ^ d 2 ( dw A (x)\ 2 d 2 W A (x) 

, 2 +U A1 (x 



dx 2 V dx I dx 1 

d 2 



dx z 

2 21 16 96 16 480 9 



U A1 {x) = x + — - + — — — - - a , + 



16x 2 3 + 4x 2 (3 + 4x 2 ) 2 15 + 4x 2 (15 + 4x 2 ) 2 2 



d 2 

dx 2 "°^"-'4> 



U(x;l)-2— log<fc{x 



This shows clearly that the non-generic two-indexed polynomial system (A) generated by a 
non-virtual state solution 0i H (x; |) and a virtual state solution 2 (x; |) is not shape invari- 
ant [10J, in contrast with the generic multi-indexed polynomial systems [2J. The standard 
Darboux-Crum transformation removes the groundstate "created" by the non-virtual state 
solution <f) l i l (x; |) and the generic one-indexed polynomial system generated by the virtual 
state solution (fi\(x; |) is obtained. 

3.3.2 Global solutions (B) 

The global solutions are 

e -v/2 n i/s 



n u <>»U) = £„<>»(.,■). n u = ~ + U B {x) ) £ n = 4n, (3.43) 



l 3 + l v) 2 C n(v), » = ().\ (3.11) 
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Ill which £ B (rj) (n > 0) is a degree n + 3 polynomial in rj defined by 

£ B ( V ) = (-63 + 25277 + 240t? 2 + dA V 3 )L^ 4 \r]) 

- 4n(3 + A7]) 2 L^ 1/4) (t]) - 4r?(3 + 4r/)(15 + Ar])d v L^ 1/A) (r]). (3.45) 

There is another member of this family corresponding to the eigenvalue £ = —12, be- 
cause of the fact that the Darboux transformation in terms of 2 n ( x ) |) is not one to one, 
(l/0aV; z), V^sPO; |)) < oo, fl3T25|) . see 0. Let us denote it tentatively by C B 3 (r]) 

C B M ^ 1. (3.46) 

Its normalisation is irrelevant. Note that n = —2, —1 are missing. Some lower members are 

297 

C B (n) = -63 + 252r] + 240t/ 2 + 64t/ 3 , £f (77) = — — + 396t? + 264r/ 2 - 64r/ 3 - 64r? 4 , 

£f (r? ) = + IZ^ + 225^ - 30077 3 - 120r? 4 + 32 V 5 , 

d , , 13167 13167 627 2 „,,„«, 4 „ _ . 32 fi 

£f0?) = -^28" + -XT' 1 + " ~ 1UV + 11271 ~ Y V ' ' ' ' ' 



They are also discussed in §6.2.6 of [9]. They satisfy the orthogonality relation 

-£%(ri)£ B (ri)dri = h B 5 nm , n, m = -3, 0, 1, . . . (3.47) 



(3 + 4^ 

h B = — ' ; -y ~ v ~ ' 4/ , n ^ -3. (3.48) 



Bd _e f 4(n + 3)(4n + 7)r(n + |) 



77 ! 



which is obtained by rewriting ( |2.10p with £l u (j) = —12 and £{(§) = —7. 
By rewriting the Schrodinger equation for 4> B {x) f!3.43j) in terms of 

<j) B {x) = e-" /2 7/l?/(7/), 

we obtain a second order differential equation with regular singularity at 77 = and 77 = — | 
and an irregular singularity at x = 00: 

„ ,3 s , 45 - 24r7+ I677 2 

vv +{- A -v)y + — (3 + 4 )2 v + ny = o, n = -3,0,1,...,. (3.49) 

For each n we have a global solution with a p = —2 apparent singularity: 

^) = 7^p n = -3,0,l,...,. 
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The differential equations for the polynomials {C B (r])} read 

477(3 + 4ri)C% (77)" - (-9 + 6477 + 16r/ 2 )£^(r/)' 

+ 4(3 + n)(3 + 4ri)C*(7i) = 0, n = -3, 0, 1, . . . , . (3.50) 

The groundstate wavefunction 4> B 3 (x) (I3.44j) with £ B 3 (13461) and the energy £ = — 12 
is related to the potential Ub(x) (13. 13p through the prepotential Wb{x) (Riccati form of 
Schrodinger equation): 

In other words, the Hamiltonian T-L B (13431) is written in a factorised form: 

n B = A B] A B - 12 A B - d dwB ^ jb\ = _^_ dw B (x) 



dx dx dx dx 

It is easy to verify that the standard Darboux-Crum transformation [5] of deleting the 
groundstate will lead to the system of one-indexed (or the exceptional) polynomial system 
generated in terms of the virtual state solution 4>\(x; |): 

„,R1 def „ r ,m d 2 ( dWB(x)\ 2 d 2 WB(x) 

H m = A B A m - 12 = -— + — f^-L - 12 

dx z \ dx J dx 1 

d 2 

+ U B i(x), 



dx 2 

2 _5_ _L6 96 _ 7 

Ubi(x) - x + + 3 + Ax2 (3 + A%2) 2 2 

= U{x-\)-2^\og]>\{x-\). 



3.3.3 Global solutions (D) 

The global solutions are 



U u <\> u n [x)=£ n <\> u n {x), n u = -^ + U D (x), £ n = 4n, (3.51) 

dx z 



- V /2 1/8 

^ (a?) = 4(-3+\)^ " (l?) ' ^ = ' 1 "--' ( 3 - 52 ) 
in which C B {rj) (n > 0) is a degree n + 3 polynomial in 77 defined by 

^(^)= f (3 + 4r / )(63 + 16r / 2 )4 7 / 4 )(r / ) 



16 



- 16m?(-3 + Ar]) 2 L^ /4) (rj) - 36r/(-3 + 4r/)(l + Arj) d v L { ^ (77). (3.53) 

There is another member of this family belonging to the energy £ = —8, to be denoted 
tentatively by £® 2 (v) 



£? 2 (77) = 1 + 477. (3-54) 



Some lower members are 



2079 

£D(rj) = (3 + 4r/)(63 + 16r/ 2 ), C?{rj) = —— + 792r/ 2 - 384r/ 3 + 192?? 4 , 

C°(ri) = - + 3465r/ 2 - 2940r/ 3 + 1512r/ 4 - 224r/ 5 , 

32 8 

„„, x 197505 65835 78375 2 , „ 4 = 352 , 

Cffr) = -^2§ 1^ + "l - 77 " 1045077 + 6270r/ " 1584r ? + ' " ' " 

Note that n = — 1 is missing. The lowest member CP_ 2 (I3.54p with the energy £ = —8 is 
related to the potential Ue>(x) ( I3.17P through the prepotential (Riccati form of Schrodinger 
equation): 

„w D (x) def e~" 2 / 2 x4(l + 4x 2 ) / dw D (x) \ 2 d 2 w D (x) _ TT 

~ (_ 3 + 4x 2 ) 2 " V dx ) + dx 2 8 -^W- 

By rewriting the Schrodinger equation for <p^(x) (I3.5ip in terms of 

we obtain a second order differential equation with regular singularity at r\ — and 77 = +| 
and an irregular singularity at x = 00: 

,. .3 . . 27 + 72r/ - 16r/ 2 
VV +{^-V)y (_ 3 + 4r/ )2 — 2/ + n 2/ = > n = -2,0,1,...,. (3.55) 

For each n we have a global solution with a p = —2 apparent singularity: 



-3 + 4r/)^ 



y»fo) = 7 o :^ 2 > n = -2,0,l,...,. 



The differential equations for the polynomials {£^(77)} read 
4r/(-3 + 4r/)C(r/)" - (1 + 4r/)(9 + 4r/)C(r/)' 



+ 4(3 + 1277 + 7i(-3 + 477))££(77) = 0, n = -2, 0, 1, . . . , . (3.56) 
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3.3.4 Global solutions (E) 

The global solutions are 



U E ^{x)=£ n ^{x) 1 n E = -^ + U E (x), £ n = 4n, (3.57) 



15 

€i*)= {6 + r]n u + r] f n(v), n = 0,l,..., (3.58) 

in which C E {rj) is a degree n + 3 polynomial in rj defined by 

C E (v) = -24(840 + 280r/ + 30r/ 2 + rf)L^\r]) 

- 4n(6 + 77) 2 (14 + v)L { J ] (v) ~ 1677(6 + rj)(12 + rj)d v L^(ri). (3.59) 

This belongs to the generic case and the Darboux transformation is one to one. The ground- 
state corresponds to C E . Some lower members are 

C E {rj) = -24(840 + 280r/ + 30r/ 2 + r/ 3 ), (3.60) 

Cf(ri) = 28(-6336 - 1320r/ + 44r/ 2 + 22r/ 3 + r/ 4 ), 

£f( v ) = -16(54432 + 4536r/ - 1944r/ 2 - 180r/ 3 + 12r/ 4 + r/ 5 ), 

C E (r]) = 6(-524160 + 21840r/ + 26520r/ 2 - 156r/ 3 - 338r/ 4 + r/ 6 ), .... 

They satisfy the orthogonality relation 

p—r\ 7 

£n(v)£l(v)dV = h n$nmi 71,771 = 0,1,... (3.61) 



o (6 + r/) 4 (14 + r/) 



/if = f 16(n + 10)(ri + 6)(n + l) 7 , (3.62) 



which is obtained by rewriting f l2TT0l) with 8\Q% ) = -40 and £f{^ ) = -24. Here (a)„ = f 
Ofc=i( a + A; — 1) is the shifted factorial (the Pochhammer symbol). 
By rewriting the Schrodinger equation for <p E (x) (13.571) in terms of 

0f (x) = e-"/ 2 r/T y (r / ), 

we obtain a second order differential equation with regular singularity at rj = and r\ = 
—6,-14 and an irregular singularity at x = oo: 

,, , , . 1008 + 1277 - 16r/ 2 - 77 3 , 

VV" + (8 - V)V' + 4 + / ?7)2 y + ny = 0, n = 0, 1, . . . , . (3.63) 
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For each n we have a global solution with a p = — 2 and a p = — 1 apparent singularities: 

^) = (6+ gi4+,,r " =0 ' 1 ''"'' 

The differential equations for the polynomials {C F (77)} read 

77(6 + 77)(14 + v)£n(v)" ~ ("672 - 877 + I877 2 + 77 3 )£f (77)' 

+ (-224+ 18r7 + 377 2 + n(6 + r7)(14 + 77))£f(r7) = 0, n = 0, 1, . . . , . (3.64) 

The groundstate wavefunction <Pq(x) (13. 58j) with C F (13.601) is related to the potential 
U E (x) (I3.19P through the prepotential w E (x): 

w„M def e- x2 / 2 xT(840 + 280x 2 + 30x 4 + x 6 ) { dw E (x)V d 2 w E (x) TT . . 
v J = r . ^ 1 =>■ . — H , „ = UfAx). 



(6 + x 2 ) 2 (14 + x 2 ) V dx J dx 2 E 

By construction, the system is shape-invariant: 

dw E (x)\ 2 <Pw E (x) _^.m_ 2 ^ logW[te m^n (:r; m ](:c)+4 



(•/./■ J dx 2 2 dx 2 2 2 

9 255 48 
x 2 + — - 14 



2 



4a; 2 (6 + x 2 ) 2 6 + x 

640(630 + 175a; 2 + 12x 4 ) 4(-160 + 3a; 2 ) 



(840 + 280a; 2 + 30a; 4 + a; 6 ) 2 (840 + 280a; 2 + 30x 4 + a; 6 ) ' 

As shown clearly above, the singularity at x 2 = —6 of the new potential corresponds to the 
characteristic exponents p = —1,2. 

3.3.5 Global solutions (F) 

The global solutions are singular in (0, 00): 

n F ^(x)=£ n ^(x), H F = -^ + U F (x), £ n = An, (3.65) 

aar 

e~ v ' 2 ri~ 4 

^0*0 = ( _ 6 + r/)2( _ 14 + r?) £ n(^) ; n = 0,l,..., (3.66) 

in which C F {jf) is a degree n + 3 polynomial in 77 defined by 

£n(v) = -4(-9(-280 + U0t] - 22r] 2 + t? 3 ) + n(-U + r/)(-6 + 77)) L^ 7) (77) 

- 16t7(-6 + 77)(-12 + 77) d v L { - 7) (77). (3.67) 
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They are not square integrable and thus they are not eigenfunctions. Some lower members 

are 

(rj) = 36(-280 + 140r7 - 22r/ 2 + r/ 3 ), (3.68) 
Cf (77) = -32(-1512 + 504ry + Ylif - 16?7 3 + r/ 4 ), 
C% (77) = 14(-6480 + IO8O77 + 396r/ 2 - A2rf - Ylrf + r/ 5 ), 
£f ( v ) = -4(-20160 + 1632r7 2 + 152r7 3 - 52r? 4 - 10r/ 5 + r? 6 ), .... 

By rewriting the Schrodinger equation for 4>^(x) (I3.65P in terms of 

13 

we obtain a second order differential equation with regular singularity at 77 = and 77 = 6, 14 

and an irregular singularity at x = 00: 

„ , s , 1008 - 12r? - 16r? 2 + 77 3 , 
W" ~ (6 + V)V' ~ 4 ( _ 6 + ??) 2 ( _ 14 V + = °> n = 0,l,...,. (3.69) 

For each n we have a singular global solution with a p = —2 and a p = — 1 apparent 
singularities: 

£ F (n) 

yn(v) = 7 ^TT 1 7, 71 = 0,1,...,. 

ynyn (-6 + r/) 2 (-14 + 77)' 
The differential equations for the polynomials {jC-niv)} read 

?j(-6 + 77) (-14 + v)£n(v)" - (504 - 1047/ - 8r/ 2 + r/ 3 )£^(r/)' 

+ (-252 - 8r/ + 3r/ 2 + n(-6 + 77)(-14 + r/))£^(77) = 0, n = 0,l,...,. (3.70) 

The wavefunction 4>o( x ) A3.66[) with £^ (I3.68|) is related to the potential Up(x) (I3.2ip 
through the prepotential wp(x): 

Wf{x) def e-" 2 / 2 x~f (-280 + 140x 2 - 22x 4 + x 6 ) f dw F {x) \ 2 d 2 w F (x) = 

(-6 + x 2 ) 2 (-14 + x 2 ) ^ V ^ / ^ 2 ~ ^ 

By construction, the system is shape- invariant. The equations and the wavefunctions of (E) 
and (F) systems have very similar forms. 

3.4 Three-index case 

Among various possibilities of the three index cases with lower degrees, the following has 
been noted to have a cubic zero: 

-x 2 /2 

W[4>[(x; f ), 4>?(x; f ), ~${x- f )](x) = -4^^(30 + x 2 ) 3 (390 + 39x 2 + x 4 ). (3.71) 

x 2 
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This belongs to the generic case. The method of construction of the global solutions is about 
the same as that of the case (E). But the explicit formulas become quite lengthy and they 
will not be reported here. 

At present we do not have a systematic means to locate seed solutions with high multi- 
plicity of zeros, we stop our investigation here. 

4 Summary and Comments 

Several families of infinitely many global solutions are presented for certain second order 
differential equations ( I3.4(jp . ( I3.49p . (I3.55p . (I3.63P and ( I3.69P having two or more finite regular 
singularities and one irregular singularity at infinity. The characteristic exponent of the 
global solutions at each finite singularity is —2, —1 or 0. They are obtained from the non- 
generic as well as the generic two-indexed Laguerre polynomials by choosing the parameter 
g in such a way that the seed solutions have a triple zero. In two families (D) and (F), the 
extra finite singularity is located inside the domain of the radial oscillator potential (0, oo). 
In the other three families, the extra singularity is outside of (0, oo) and thus they form 
orthogonal polynomials over (0, oo) with the weight functions (A) Wa(ji) = e~ r, r] 1 / 4 j (3+4^) 4 , 
(B) W B (v) = e~V 1/4 /(3 + 4r])\ (E) W E (v) = e~V/{(6 + r/) 4 (14 + r/) 2 }. Since they are 
obtained within the framework of multi-indexed orthogonal polynomials, their degrees have 
"gaps". Those polynomials belonging to (A) and (B) were derived in a different way in [9]. 

By construction, these families of two-indexed Laguerre polynomials are also the main 
part of eigenf unctions of exactly solvable quantum mechanics, belonging to the class called 
deformed radial oscillator potentials. In contrast with the generic multi-indexed polynomials 
(E) and (G), the systems obtained from the non-generic multi- indexed polynomials (A) and 
(B) are not shape invariant. By the standard Darboux-Crum transformations, these non- 
generic two-indexed polynomial systems are mapped to the generic one-indexed polynomial 
systems. 

By using the same method, we tried in vain to find other global solutions with higher 
degrees (p = —3, —4) of apparent singularities within the multi-indexed Laguerre polynomi- 
als. It would be interesting to explore the Jacobi polynomials counterparts of those results 
established in the present paper [11]. We wonder if the ideas of the present paper could be 
generalised to a certain kind of difference equations in discrete quantum mechanics [T2]-[18j. 
The multi-indexed orthogonal polynomials are also established in the framework of difference 
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Schrodinger equations, the multi-indexed (g-)Racah, Wilson and Askey- Wilson polynomials 

mm- 
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